Mon. Not. R. Astron. Soc. 000. HHT7lf 20061 Printed 5 February 2008 (MN WF$t style file v2.2) 



Impact of a non-Gaussian density field on 
Sunyaev-Zeldovich observables 

Sharon Sadeh 1 *, Yoel Rephaeli 1 ' 2 and Joseph Silk 3 

- - - 1 School of Physics and Astronomy, Tel Aviv University, Tel Aviv, 69978, Israel 

Center for Astrophysics and Space Sciences, University of California, San Diego, La Jolla, CA 92093-0424 

(N 



Of 



O 



X 



5 February 2008 



ABSTRACT 

The main statistical properties of the Sunyaev-Zeldovich (S-Z) effect - the power spec- 
. trum, cluster number counts, and angular correlation function - are calculated and 

compared within the framework of two density fields which differ in their predictions 
of the cluster mass function at high redshifts. We do so for the usual Press & Schechter 
mass function, which is derived on the basis of a Gaussian density fluctuation field, and 
£^ . for a mass function based on a x 2 distributed density field. These three S-Z observables 

' are found to be very significantly dependent on the choice of the mass function. The 

f^i . different predictions of the Gaussian and non-Gaussian density fields are probed in 

detail by investigating the behaviour of the three S-Z observables in terms of cluster 
mass and redshift. The formation time distribution of clusters is also demonstrated to 
be sensitive to the underlying mass function. A semi-quantitative assessment is given 
of its impact on the concentration parameter and the temperature of intracluster gas. 
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C3 ■ 1 INTRODUCTION 

It has recently been suggested that several observational results may be hinting at a significant non-Gaussian component in the 
primordial density fluctuation field (Mathis, Diego, & Silk 2004, hereafter MDS). Among these are the detection of structures 
with high velocity dispersions at redshifts z = 4.1 (Miley et al. 2004) and z = 2.1 (Kurk et al. 2004), the apparent slow 
evolution of the X-ray cluster luminosity function indicated by objects taken from the Vikhlinin et al. (1998) and Mullis et al. 
(2003) cluster catalogs, and the enhancement of the CMB power spectrum at high multipoles that was measured by the CBI 
(Mason et al. 2003) and ACBAR (Kuo et al. 2004) experiments. MDS pointed out that scale-dependent non-Gaussianity may 
be probed through its detectable manifestations on the large mass scales of clusters of galaxies, and that this is possible without 
necessarily violating constraints on non-Gaussianity such as those obtained from WMAP all-sky CMB maps (Komatsu et al. 
2003). It is argued by MDS that while it is hard (although possible) to reconcile the observational results with predictions 
made on the basis of a Gaussian field, these arise more naturally with either positively-or negatively-skewed density fields. 

In this paper we extend the work of MDS, focusing on the impact of a positively-skewed non-Gaussian field on the 
principal statistical properties of the Sunyaev-Zeldovich (S-Z) effect. We compare in detail the predicted S-Z power spectrum, 
cluster number counts, and the angular correlation function for Gaussian and non-Gaussian models. 

In section 2 we specify the various models and scalings used in order to carry out the calculations; specifically, a complete 
account of the mass functions and the relevant scalings for the quantification of cluster-induced Comptonization are provided. 
Results of the calculations relating to the S-Z power spectrum and number counts are detailed in section 3. In section 4 we 
formulate the method of calculation of the 2-point angular correlation function of S-Z clusters and present the corresponding 
results. The calculation of cluster formation times within the framework of non-Gaussian fields requires a slight modification 
of the corresponding Gaussian case. This is explained in section 5. Results of the formation time computations in both the 
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Gaussian and non-Gaussian models follow in the same section, followed by a brief discussion and a summary of our main 
conclusions in section 6. 



2 THE MODEL 

Calculations of statistical properties of the evolving cluster population involve the basic cosmological and large-scale quantities 
and parameters, as well as essential properties of intracluster (IC) gas. Detailed description of the calculation of integrated 
S-Z properties can be found in the literature. Therefore, our brief treatment here highlights only those aspects that pertain 
to the basic properties of the density fields explored in this work. The multiplicity function in the Press & Schechter (1974) 
formalism is 

n(M,z) = -F M £-^ d M, (1) 

where [i = S c /a and pt is the background density. The critical overdensity for spherical collapse is assumed to be constant, 
S c = 1.69, owing to its weak dependence on redshift. The density field variance, smoothed over a top-hat window function of 
radius R, is calculated through 

p oo 

a 2 (R) = / P(k)W 2 (kR)k 2 dk, (2) 
Jo 

where P(k) = Ak n T 2 (k) . The mass variance evolves with redshift according to 
where 

"»W = {1 + z) $lTto)?-nmY " A(2) = 1 - Qmiz) > (4) 

and (Carroll, Press & Turner 1992) 

9[Qm{z)] = [n m ( z yn - n A + (i + n m ( z )/2)(i + n A /70)] ' (5) 

The function F(fi) in equation Q assumes the form 

FM = Jl e -4jL ) (6 ) 

V 7T PR 

for a Gaussian density field, and the form 

for a density fluctuation field which is distributed according to the Xm model (Koyama, Soda & Taruya 1999) . The parameter 
m describes the number of CDM fields added in quadratures, so as to yield the requested distribution. 

In most of our calculations we use m = 1, which constitutes the largest possible deviation from a normally distributed 
probability density function (PDF) for this model. Increasing the number of random density fields obviously results in a PDF 
that approaches a Gaussian, in accordance with the central limit theorem. Some of the calculations were done also for the 
case m = 2 in order to illustrate the impact of this parameter on the high-mass tail of the PDF. The respective mass functions 
are then obtained by multiplying either expression by a common factor, as explicitly stated in equation Q. 

For the Gaussian case an adiabatic CDM transfer function is employed, 

ln(l + 2.34g) 

[J. 1- o.o»j -r ^iu.i(j; T (O.wq) -|- yu.i it, 
whereas an isocurvature CDM transfer function is chosen for the Xm model 



T{k) C DM = l 2 34g q> [1 + 3-89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ]- 1/4 , (8) 



T(k)cDM,isoc = (5.6g) 



(40g) 2 . ws 1 " 

1 + 215g + (16g) 2 (l + O.Sg)- 1 1 q) 



0) 



with q = k/(fl m h 2 ) Mpc -1 . Both transfer functions were taken from Bardeen et al. (1986). The calculations were carried out 
for a ACDM model, with Ha = 0.7, fi m = 0.3, h — 0.7. In the Gaussian case the spectral index was taken to be n = 1 and the 
normalization as = 0.9. In the non-Gaussian case n = —1.8, and from the requirement that the present cluster abundance is 
the same as calculated in the Gaussian model, we obtain as = 0.73. 

The S-Z angular power spectrum is calculated using the basic expression 
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Ce= f r 2 ^- [ n(M,z)G e (M,z)dMdz, (10) 
Jz dz Jm 

where r is the co-moving radial distance, and Gt is obtained from the angular Fourier transform of the profile of the temperature 
change induced by Comptonization of the CMB at an angular distance 9 from the centre of a cluster, AT (6). The function 
Gi, which is proportional to AT(6) 2 , is fully specified in, e.g., Molnar & Birkinshaw (2000). In the limit of non-relativistic 
electron velocities (more on this in section 3.2) the relative temperature change assumes the simple form 

^ W =[xcoth(|)-4]j / W, (11) 
where, for an isothermal cluster with a /3 density profile, the Comptonization parameter is 



(<9) = 2 kB ° T no( - M ' z ) r °( M ' Z M M ' z ) tan -i 



1 + (6/dc) 2 



(12) 



Here no, To, r c , and 6 C denote the central electron density, temperature, core radius, and the angle subtended by the core, 
respectively. The virial radius is taken throughout to be p — 10 times r c . 

Of crucial importance are the scalings of cluster properties with mass and redshift; the temperature is assumed to scale 

as 

where a temperature of 8.5 keV is ascribed to a cluster with mass 10 15 /i _1 Mq at present. The parameter a is usually taken to 
be 2/3, in accordance with theoretical predictions based on hydrostatic equilibrium. We allow a variance of 10% in its value, 
as to possibly address non-gravitational effects which may result in a deviation from the simplistic assumption of hydrostatic 
equilibrium, and also to reflect the observed variance in this parameter. An additional uncertainty relating to this scaling is 
represented in the parameter ip, which describes the temporal evolution of cluster temperatures, and is taken to be either 
zero or unity, constituting two limiting cases for which there is either no evolution, or strong evolution, respectively. Little 
is known on the evolution of the X-ray temperature function due to the lack of X-ray observations of high-redshift clusters, 
although some interesting insight is gained from hydrodynamical simulations (Norman 2005) . 
The electron density is parameterized as 

n e {M 1 z)~n -L{l + z)\ (14) 

with the scaling of the gas mass fraction to ~ 10% (e.g. , Carlstrom, Holder & Reese 2002). Since we do not yet know the 
redshift dependence of /, we assume / = 0.1 to be roughly valid throughout the redshift considered here. Finally, the core 
radius is calibrated according to the simple relation 

M \ 1/3 1 

10^-imJ — MPC - (15) 

In accord with observational results, we adopt a variance of 20% in the value of the core radius. In order for the mass to 
remain constant when different scalings are used for the core radius, a corresponding change is affected in the central gas 
density. 

The number of clusters with S-Z flux (change) above AF„ is (e.g. Colafrancesco et al. 1997) 

iV(>AF„) = [ r 2 ^-dz [ _ n(M,z)dM. (16) 
J dz Jaf„ 

For a cluster with mass M at redshift z 

where 'je and 7 denote line of sight (los) directions through the cluster centre, and relative to this central los, respectively, 

y = 2^^n {M, z)T {M, z)r c (M, z), (18) 
m e cr 

and the spectral dependence (of the thermal component) of the effect is given (in the nonrelativistic limit; more on this in 
section 3.2) by 

= ( e »-i)2 b coth - 4 ' ( 19 > 
where x = hv/kT. The profile of the effect is given in 
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y(\ie\,M,z) 



Vi + {o/6c 



tan 



i - {6/ P 8 c y 



1 + (6/0c 



(20) 



a B ) denotes the angular 



(21) 



which is the los integral along a direction that forms an angle 6 with the cluster centre. -R s (|7 — ye\ 
response of a detector whose beam size is given in terms of ob- Finally, 

_ jAF v E{u)du 
fE(u)du ' 

is the flux weighted over the detector spectral response E(v). 

Only clusters for which the S-Z flux exceeds the PLANCK/HFI detection limit of 30 mjy are considered here. This is 
manifested in the lower limit of the mass function integral, equation 11611 . The results for the S-Z power spectrum and number 
counts described below refer to these two models: (a) IC gas temperature evolves with time according to ip = 1, and (b) no 
temperature evolution, ip = 0. Power spectra were calculated at v = 353 GHz, with a beam size is 7.1'. Number counts were 
calculated also at v = 143 GHz and v = 545 GHz. 



3 RESULTS 

3.1 Power spectrum 

Several qualitative assessments can be pointed out before we describe specific results. The longer tail of the non-Gaussian 
distribution function (with respect to that of the Gaussian) at the high mass end translates to higher cluster abundance 
in general, and high mass clusters in particular. This holds in the entire redshift range (with the exception of z = 0, since 
the mass function in the non-Gaussian case was normalized such as to yield the same cluster abundance as in the Gaussian 
case). This tendency stems from the fact that a positively-skewed density fluctuation field exhibits earlier cluster formation 
by enhancing the probability that the average density within a spherical region of any size exceeds the critical density for 
collapse. Consequently, the S-Z power spectrum is expected to attain higher levels in the non-Gaussian case. Owing to the 
inferred denser population of clusters at high redshifts, the magnitude at the peak should shift to higher multipoles, reflecting 
the higher density of distant clusters with smaller angular sizes. Naturally, this effect would be drastically reduced in a non- 
evolving temperature scenario due to the inferred lower temperatures associated with high redshift clusters when the (1 + z) 
scaling is factored out. 

The S-Z power spectrum calculated for cases (a) and (b) is plotted in Figs. through Qj] (with the exception of Fig. [3J|. 
In all plots the upper and lower panels correspond to cases (a) and (b), respectively, whereas the left- and right-hand panels 
pertain to the Gaussian and Xi models, respectively. Partial contributions to the total power spectrum from redshift intervals 
of Az — 1 are shown in Fig. Q Fig. depicts partial contributions from (logarithmic) mass intervals of log AM = 1 to the 
total power spectrum. Power variation corresponding to a 10% variance in the parameter a, the scaling of cluster temperature 
with mass, is illustrated in Fig. 0] Finally, the impact of introducing a variance of 20% in the cluster core scaling on the 
magnitude of the S-Z power spectrum is shown in Fig.|K| 

Referring first to the magnitude of the total power spectrum in Fig. corresponding to the redshift interval < z < 6, 
in both cases (a) and model (b) higher levels of power are clearly evident in the non-Gaussian model. Quantitatively, the 
maximum power levels attained in case (a) are ~ 3 • 10~ 12 and ~ 6 ■ 10 -11 for the Gaussian and non-Gaussian models, 
respectively. Note that both peak at multipoles higher than t > 10, 000; in the Gaussian model this is due to the high gas 
densities of distant clusters implied by the constancy of the gas fraction, whereas in the non-Gaussian model the combination 
of high gas densities and the long tail at the high-mass cluster end is responsible for this result. Compared to case (a), distant 
clusters are cooler in case (b) by virtue of the redshift independence of their temperatures, resulting in lower power levels 
(~ 8 -13 and 5 ■ 10 -12 ) in the Gaussian and x 2 models, respectively. In the Gaussian model the power peaks at I ~ 6000, 
reflecting the lower contribution from cooler, more distant clusters. In the non-Gaussian model, a relatively high abundance 
of massive clusters at high redshifts (i.e, small angular scales) leads to sustained high power levels on these scales. It is also 
worth recalling that regardless of whether or not the temperature changes with redshift, it still increases with increasing mass, 
so even when tp = a denser population of hot clusters is produced in the long tail of the x 2 PDF. 

The contribution of clusters lying at relatively high redshifts to the overall power is clearly seen in the \ 2 model. In fact, 
in case (a) the contribution of clusters lying in the redshift range < z < 1 peaks at t ~ 3000. However, the total power 
continues to rise owing to the non-negligible power generated by clusters situated at z > 1, particularly so from redshift range 
1 < z < 2. In fact, at £ ~ 8000 this contribution is surpassed by the next redshift range, 2 < z < 3. In the Gaussian model the 
contribution from this range to the overall power is decidedly negligible. The same effect, although less pronounced, is evident 
in case (b); here distant clusters are cooler [than in the corresponding case (a)] by virtue of the ip = scaling. Consequently, 
the partial contributions from redshift ranges 1 < z < 3 is still discernible, but to a lesser degree than in case (a). 

Fig. [5] illustrates the distribution of S-Z power in four mass ranges. Most striking is the relative contribution to power 
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Figure 1. Angular power spectrum of the S-Z effect. Shown are levels of power contributed by clusters lying at redshift range of < z < 6, 
and at 6 redshift intervals. Left and right-hand panels correspond to the Gaussian and \\ models, respectively. Upper and lower panels 
relate to cases (a) and (b), respectively 



from clusters lying in the 1O 15 /i _1 M < M < 1O 16 /i _1 M mass range. Whereas it is practically negligible at all multipoles in 
the Gaussian model, this is not the case in the non-Gaussian model; in fact, these massive clusters dominate the total power 
up to t ~ 1000 and I ~ 2000 in cases (a) and (b), respectively. At higher multipoles the 10 14 /i _1 M Q < M < 1O 15 /i _1 M 
range dominates. In the Gaussian model low-mass clusters of 10 13 M < M < 1O 14 /i _1 M prevail only at the highest 
multipoles. These clusters are obviously the most abundant objects predicted by a mass function based on a primordial 
density fluctuation field of Gaussian nature, and - in light of their relative small sizes (and low temperatures) - thus must 
lie at a relative proximity to the observer; more distant clusters of similar masses would not produce sufficiently strong S-Z 
signal to be detected. Otherwise, in the Gaussian case, the lion's share of power originates in clusters in the mass range 
lO^/i^Mg < M < W 15 h- 1 M Q . 

Fig.0and|5|explicitly demonstrate the significantly higher contribution of high-redshift - and therefore hotter, in particular 
in case (a) - clusters in the x? model. This is manifested in higher power magnitude and at an intensification at higher 
multipoles corresponding to distant, low angular-size clusters. To further contrast the predictions of the two models, the ratio 
of total level of the S-Z power in the xi model to that in the Gaussian model for cases (a) and (b) is shown in Fig. El Fot 
I = 10 3 — 10 4 , this ratio increases from 7 to 20 in case (a), and from ~ 4 to ~ 7 in case (b). 

The effect of changing a, the scaling of the cluster temperature with mass, is shown in Fig.|I] Since T oc (M / Mis) a , two 
different behaviours are expected for M < Mis and M > Mis. For M < Mis, increasing a with respect to the conventional 
value of 2/3 results in cooler clusters, whereas if M > Mis, higher temperatures will be attained. Lowering a results in the 
opposite effect, i.e., warmer low-mass clusters and cooler high-mass clusters. These changes affect the power spectrum for the 
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Figure 2. The same as in the previous figure, but here power spectrum levels are calculated for 4 mass intervals, spanning the range 
lO 13 ft- 1 M < M < lO^/i-^g. 
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Figure 3. A comparison of S-Z power spectrum levels in the X\ an d Gaussian models for cases (a) and (b). 
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Figure 4. The impact of the temperature-mass relation on the S-Z power spectrum, T oc (M/Misp. Shown are 3 cases: £ = 0.7(3), £ = 
0.6(6), § = 0.6. The arrangement of the panels is the same as in the previous figures. 



Gaussian and non-Gaussian models in slightly different ways; we have already seen that in the Gaussian model the largest 
contribution to the power spectrum comes from clusters with masses 10 14 /i _:L M Q < M < 10 15 /i _1 M . That is, an outcome of 
increasing a is that most of the contributing clusters undergo cooling, and the level of power decreases. Obviously, lowering 
a has the opposite effect. This is clearly seen in the upper- and lower- left panels of Fig 2] It can also be seen (particularly 
so in the lower-left panel) that the curves peak at decreasing multipoles with increasing a. This has to do with the fact that 
the contribution of these relatively low-mass clusters to the overall power levels decreases with increasing a, so that power is 
suppressed from high multipoles, resulting in a leftward shift of the peaks. In the Xl model the results are somewhat different. 
Here the contribution to the power from high-mass clusters (10 15 /i _1 M Q < M < lO 16 ft _1 M ) dominates at £ < 1000 (case 
a) and £ < 2000 (case b), whereas the corresponding contribution at higher multipoles originates in clusters lying within the 
mass range 10 14 /i _1 Mq < M < 10 15 /i _1 Mq. In other words, an increment of power with increasing a is expected at lower 
multipoles, and a decrease of power with increasing a at higher multipoles. This is shown in the upper- and lower-right panels 
of Fig. m 

We also investigate the influence of r c on the power spectrum levels. The results are illustrated in Fig. and are rather 
obvious. First, lowering or increasing the core radius changes the effective cluster size, and consequently, the multipole at 
which the power peaks. This can be easily seen in the lower-left panel of Fig. [21 In the remaining panels it is somewhat harder 
to notice the effect, owing to the fact that the peaks fall beyond £ — 10000. Also, in order to keep the total gas mass constant, 
lower and higher core sizes necessitate higher and lower central gas densities, respectively. This results in higher and lower 
power levels for lower and higher core-sized clusters, respectively, as is apparent from the fact that the power spectrum scales 
as the central gas density squared (equations 1101 and I12H . 
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Figure 5. The impact of cluster core scaling on the S-Z power spectrum. Each of the three values considered (r c 
0.1, 0.125, 0.15 Mpch -1 ) denote the core radius of a cluster with mass 10 15 Mq h~ 1 at redshift 2 = 0. 



It is useful to mention some recent studies of S-Z power spectrum calculations and compare with our results. Case (a) 
of the Gaussian model is the most adequate for comparison, since past work has focused on Gaussian density fields, and 
whether analytically or inferred from hydrodynamical simulations, the IC gas temperature has been assumed (or observed) 
to increase with redshift (at least to a certain degree), in contrast with the constant temperature scenario, case (b). Da Silva 
et al. (2001) employed hydrodynamical simulations to find power levels of 5 — 9 ■ 10 12 peaking at £ ~ 7000, depending on 
whether or not non-adiabatic processes such as preheating and radiative cooling were included. Their selection of cosmological 
parameters is essentially identical to ours, but the calculations were carried out in the Rayleigh- Jeans region, so a comparison 
between their results and ours requires an introduction of a multiplicative factor of 4 to our power levels. Consequently, the 
magnitude of the total power for case (a) in the Gaussian model peaks at multipole I ~ 10, 000, at a value of ~ 1.2 ■ 10 -11 . 
Bond et al. (2005) have quantified the S-Z power spectrum using both analytic calculations and results of hydrodynamical 
simulations. The simulations yield peak power levels of 1 ~ 4 ■ 10 -11 at £ ~ 10,000. Their analytic calculation yields similar 
power levels, but at lower peak multipoles, £ ~ 5000 — 6000. Note that they calculated S-Z power at the CBI frequency, 
v = 31 GHz. Rescaling our results to this frequency gives a power level of ~ 1.1 ■ 10~ n . Refregier & Teyssier (2002) have 
also used hydrodynamical simulations and an analytic approach to calculate the S-Z power spectrum. They found that both 
methods yield approximately the same power levels (peaking at £ ~ 8000 with magnitude > 10 -11 ), particularly so when a 
cluster mass range 5 ■ 10 10 < M < 8 ■ 10 14 h -1 Mq is employed in the analytical calculation. The calculations were carried 
out in the R-J region. Springel, White & Hernquist (2001) employed hydrodynamical simulations, neglecting non-adiabatic 
processes, to infer power levels of 2 — 4- 10 -11 at the peak, £ ~ 10, 000. In fact, the latter paper presents a detailed comparison 
of S-Z power spectra published in the literature showing a large scatter (spanning over an order of magnitude), which can 
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Figure 6. Total angular power spectrum of the S-Z effect as calculated above, together with results from the CBI and ACBAR 
experiments. 

largely be attributed to uncertainties in the theoretical modeling of IC gas and the techniques implemented in hydrodynamical 
calculations. To complete this section, we combine, in a single plot, the total power spectrum levels of the S-Z effect in the 
Gaussian and Xi models, cases (a) and (b), with results from the CBI and ACBAR experiments. The analytic results were 
rescaled so as to correspond to an observation frequency of 31 GHz, corresponding to the CBI instrument. As can be seen 
in Fig. ©, it seems difficult to reconcile the Gaussian results with the data, whereas case (a) in the \i model provides a 
reasonable match to the observed power excess at £ > 2000. These results agree with those found by Mathis et al. (2004). 



3.2 Number counts 

The qualitative considerations pertaining to the power spectrum that were presented in the previous section apply also to 
S-Z number counts: Higher (cumulative) counts are expected in the non-Gaussian case, and those are likely to be distributed 
differently in redshift space with respect to the Gaussian case, reflecting the higher population of massive clusters at high 
redshifts, which satisfy the detection flux limit. Fig. |7|depicts the redshift distribution of cumulative S-Z number counts. The 
arrangement of plots according to model is the same as in Fig. Q through ©. In none of the models clusters with masses 
lower than lO 14 /i _1 M generate fluxes that exceed the adopted detection limit, 30 mjy. 

The most noticeable results displayed in the figure are: In the Gaussian model the largest contribution to the cumulative 
number counts (at z — 0) comes from clusters lying in the mass range lO 14 /i _1 M < M < 10 ft Af©. Actually, the 
cumulative number counts at z — in this mass range are higher by at least an order of magnitude than those contributed 
by clusters in the mass range 1Q 1S h^MQ < M < 1O 18 /i _1 M . Specifically, in case (a) the curves depicting the total counts 
and those contributed by clusters with masses in the range 1O 14 /i _1 M < M < 10 15 hT 1 Mq are virtually identical. In case 
(b), however, at redshifts z > 0.5 counts are dominated by clusters with higher masses, 

This can be attributed to the fact that in case (b) the clusters in the former mass range have lower S-Z fluxes due to the 
lower temperatures as compared to clusters in case (a). Clusters belonging to the latter mass range are still massive enough 
(and consequently, hot enough) to compensate for the lack of temperature evolution. Regardless of this effect, it should be 
emphasized that the counts in the high mass range are negligible in comparison with the total cumulative counts at z = 0. 
A slightly different behaviour is apparent in case (a) of the x\ model; here total counts are still dominated by clusters in 
the mass range lO 14 /i _1 M < M < 10 15 /i _1 M Q , although differences are not as pronounced as in the Gaussian case. While 
counts in the lower-mass range in the Gaussian case are higher by a factor of ~ 30 — 40 than the corresponding counts in the 
higher-mass range, here a mere factor of ~ 3 is seen. In case (b) the results are even more interesting, since the total counts 
are dominated by high-mass clusters, lying at < z < 3.5. The cumulative counts at the lowest redshifts consist of a slightly 
higher contribution from the high-mass range and are larger by a factor ~ 1.3 than the corresponding contribution from mass 
range 10 14 /i _1 M < M < 1O 15 /i _1 M . Differences in total counts are large; in the Gaussian model these are ~ 4000 and 
~ 3000 in cases (a) and (b), respectively, as compared to levels of ~ 30, 000 and ~ 10, 000 in the corresponding cases of the 
Xi model. 

Fig. |H| describes cumulative number counts for two more frequency channels, v = 143 GHz and v — 545 GHz, in 
addition to v = 353 GHz which was used above. In the non-relativistic electron velocity limit the (thermal component of 
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Figure 7. Cumulative S-Z counts as a function of redshift; shown are contributions from 4 mass intervals. Note that the contribution 
from the lowest mass interval (10 13 < M < 10 14 solar masses) vanishes in all cases considered. 



the) S-Z spectral function - which is essentially independent of the gas temperature - is g(x) (eq. 19). At these frequencies, 
g(x) = —4, 3.2, 6.7, respectively, so our results for v = 353 GHz can be easily scaled accordingly. Adopting the same value of 
the limiting flux (30m Jy), we find that values of the ratios of number counts at v = 353 GHz and v = 545 GHz are larger in 
the Gaussian than in the xi model. This can clearly be attributed to the significantly higher population of clusters with flux 
exceeding the limit in the x\ model. Note that an accurate calculation of the S-Z intensity change necessitates a relativistic 
calculation (Rephaeli 1995) which yields a more complicated expression for the temperature- dependent spectral function (for 
which there are a few analytic approximations; see, e.g., Shimon & Rephaeli 2004, Itoh & Nozawa 2004). The deviation from 
the nonrelativistic calculation can amount to few tens of percent for typical temperatures in the range 5-10 keV. Since our 
main focus here is a comparison between predictions of the two density fields for quantities that are integrated over the cluster 
population (rather than an accurate description of the effect in a given cluster), it suffices to use the much simpler function 
9{x)- 

We may summarize the results of this subsection as follows: total count levels are significantly higher in the xi model, 
roughly by factors of ~ 7.5 and ~ 3 in cases (a) and (b), respectively. The Contribution of high-mass clusters to the total 
cumulative counts are practically negligible in the Gaussian case, but important in the non-Gaussian model, particularly so 
in case (b) where the lower-mass clusters fail to generate sufficient S-Z flux by virtue of their lower temperatures with respect 
to the evolving temperature scenario in case (a). In addition to this, the counts are more broadly distributed in redshift space 
in the non-Gaussian model due to the higher abundance of clusters at high redshifts. 

As in the section describing the S-Z power spectrum, a comparison with published works on S-Z number counts is 
pertinent. Springel et al. (2001) find ~ 0.6 clusters per square degree for an observation frequency of 150 GHz and a flux 
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Figure 8. Cumulative S-Z counts as a function of rcdshift. Number counts are shown for the three specified frequency channels of the 
HFI instrument on the Planck satellite. 

limit similar to ours, 30 mjy. This translates into all-sky counts of ~ 2.5 ~ 10 4 clusters, whereas we find ~ 4 • 10 3 . Holder 
et al. (2000) report a prediction of ~ 25 clusters per square degree, or a total of over 10 6 clusters. However, it is not clear 
what was the flux limit they assumed. In addition, the slightly higher normalization they adopted, as = 1, may be partly 
responsible for the inferred high counts. Benson, Reichardt & Kamionkowski (2002) find ~ 0.2 clusters per square degree for 
PLANCK specifics, with observation frequency 143 GHz and flux limit 30 mjy. An all-sky survey is then expected to yield 
~ 8000 clusters. While our predictions are lower, it should be emphasized that similar to results for the S-Z power spectrum, 
S-Z cluster counts are equally susceptible to uncertainties in IC gas modeling and to the specifics of the experimental setting, 
such that a large scatter in the results is unavoidable. Moreover, our main concern in this study is the difference between 
the magnitudes of S-Z observables as predicted by two different mass functions, and those would certainly persist were we to 
modify the modeling of the IC gas and any other component relevant to the calculation of the effect. 

3.3 A comparison between the Xi an d xl cases 

It is useful to examine the sensitivity of the results corresponding to the Xm model to the number of fields added in quadratures, 
m. We have repeated the numerical calculations of the S-Z power spectrum and number counts (and in addition, the 2-point 
angular correlation function of S-Z clusters, which is dealt with separately in the next section) with m = 2 for case (a). Since 
the degree of skewness of the initial probability density function is reduced with increasing m, this will be manifested in 
the mass function as well, and necessitates a different normalization parameter as. For m = 2 it is found that in order to 
reproduce the observed cluster abundance at z = 0, as — 0.77, instead of the lower value of as ~ 0.73 in the m = 1 model. 
Results of the calculations are illustrated in Fig. [5] The differences are quite significant in both plots; power levels in the 
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Figure 9. A comparison between S-Z power levels (left-hand panel) and number counts (right-hand panel) for a Xm'distributed density 
field with m = 1 and m = 2 is shown. Calculations were carried out for case (a). 



m — 1 case are higher by a factor of ~ 1.5 — 1.8 across the multipole range of 100 ^ I ^ 10000. It is the highest-mass clusters 
that are most susceptible to the high-mass tail of the mass function, which is strongly dependent on the parameter m. In 
particular, this is most noticeable at high redshifts where the predictions of the mass function for the massive clusters in the 
m = 1 and m — 2 cases vary greatly. Since these distant clusters have small angular sizes, it is not surprising that the ratio 
between the two curves increases with increasing multipole, reflecting the relative higher abundance of clusters with small 
apparent sizes for m = 1. This is further established by noting that the number counts ratio for m = 1 and m = 2 are highest 
at high redshifts: While the ratios are ~ 6 and ~ 3.5 at z ~ 4 and z ~ 2, respectively, the cumulative counts for m = 1 at 
z = 0.1 are approximately ~ 1.5 times as high as those obtained for m = 2. Clearly, very massive clusters at high redshifts 
are more common in the model with m — 1 than in that with m = 2. 



4 ANGULAR CORRELATION FUNCTION OF S-Z CLUSTERS 

The angular 2-point correlation function (hereafter A2PCF) is expected to be steeply dependent on the mass function, even 
more strongly than the angular power spectrum and cluster number counts, due to the quadratic dependence on the cluster 
abundance. The A2PCF of S-Z clusters has been studied by Diaferio et al. (2003), who explored the dependence of the 
expected number of cluster pairs on the biasing relation between the cluster and the mass distribution, and Mei & Bartlett 
(2003), who studied the feasibility of alleviating the degeneracy between Q m and as by means of combining S-Z cluster counts 
with their angular correlation function. The requisite data will become available when targeted S-Z cluster surveys are carried 
out. As pointed out by Mei & Bartlett, the scientific yield from the A2PCF lies in the fact that S-Z cluster catalogs, on their 
own, may constitute an important probe of cosmological models and cluster properties, even without the detailed follow-up 
observations of individual clusters. Here we calculate the A2PCF of S-Z clusters for the Gaussian and \i mass functions and 
assess its diagnostic use to distinguish between the two models. 



4.1 Formalism 



The angular correlation function can be calculated by equating the surface density of S-Z clusters with their integrated spatial 
density: 



r 2 M^M„(M )Z ) dzdM 
dz 



[1 + w{0)\ = dfii dQ 2 



r 2 ( Zl )^j^r 2 ( Z2 )^j^ m n 2 [1 + £ (Mi ,M 2 ,zi,z 2 , R)] dz x dz 2 dM x dM 2 , 



dzi 



dz 2 



(22) 



I Z! J z 2 J Mi J M 2 

where £ (Mi, M 2 , zi, z 2 ,R) is the spatial 2-point correlation function of clusters with masses Mi and Ms, located at redshifts 
Z\ and 22, and separated by physical distance R, and rii = n(Mi,Zi). As in the calculation of the number counts, the lower 
limits of the integrals are set so as to yield S-Z fluxes that exceed the detection limit. Since the uncorrelated surface density 
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term on the left-hand side of equation 1221 is identical with the integrated spatial density of the uncorrelated term on the 
right-hand side, the angular correlation function can be evaluated as 

' dr<yZ2> m n 2 £ (Mi, Mi, zi,Z2, R) dzi dz 2 dMi dM 2 

. (23) 



w(8) = 



Z1 J z 2 J Mi J M 2 



• .rfr(zi) 2 



r (zi)- 



dz2 



r 2 (z)^-^n(M,z) dz dM 



It remains to specify the spatial correlation function of clusters, £. This can be done either by employing the observational 
correlation function of clusters (e.g. Bahcall & Soneira 1983), or by regarding it as the biased counterpart of the correlation 
function of dark matter, £dm- We choose to employ the second method, in which the cluster correlation function is commonly 
written as 



$(M u M 2 ,zi,Z2,R) = b{M 1 ,z 1 )b{M 2 ,z 2 )i dm {R,z), 



(24) 



where b is the linear bias factor. Making use of the global isotropy of the universe, we can express the correlation function as 
the Fourier-transformed power spectrum of dark matter: 

A 

2^2 



£dm(R, Z) 



e ik - R P(k) d 3 k = -^- 



k n + 2 T 2 {k) smkR dk 



kR 



(25) 



The simplest form for the bias factor is (Matarrese et al. 1997, Catelan et al. 1998) 
_5c 1_ 

lr(z) 8 C - 



b{M,z) 



1 + 



(26) 



Note that more accurate formulae for the bias factor exist, which take into account non-spherical collapse (Sheth, Mo & 
Tormen 2001) and nonlinear clustering effects (Peacock & Dodds 1996); we will limit ourselves here to the simplest case, 
since our main objective is this study is a comparative evaluation of the predictions of two mass functions (rather than their 
respective accuracy). 

To proceed, it is necessary to write down the physical distance R between the members of a cluster pair, separated by 
an angle 6. In doing so it is very convenient to take advantage of the fact that the observational spatial correlation function 
falls off rapidly with distance, £.(R) ~ R~ 18 (e.g. Bahcall & Soneira 1983); a corresponding trend can be inferred from 
equation 12511 . although in this case it is obviously different from a power law. Consequently, correlated clusters must lie at 
low angular separations and at approximately the same redshift. The separation can be written as 



R=y/[r(z 1 )-r(z 2 )] 2 + (d A e) 2 , 



(27) 



where r and d A denote the radial and angular diameter distance, respectively. Substituting z = (zi+z 2 )/2 and u = z\~z 2 << 1 
yields 



dv dv 
r(zi) - r(z 2 ) = r(z 2 +u) - r(z 2 ) » r(z 2 ) + - — (z 2 ) ■ u - r(z 2 ) » u ■ — 

az 2 dz 

and the separation simply becomes 
R = 



dr 
dz 



+ [d A (z) ■ 9} 2 



(28) 



(29) 



Substituting this in the numerator of equation 12311 . together with the spatial correlation function of equation 12411 . we have 



dr(zi) 



\z 2 ) 



dzi dz 2 
so that we can now write 



dr ^ ~r\z)(^) 2 dudz, 



(30) 



w(9) = 



z J Mi J M 2 



r 4 (z) 



dr(z) 
dz 



n(Mi,z) n(M 2 , z) bi(Mi,z) b 2 (M 2 , z) dM! dM 2 dz 



[^■u] 2 + [d A (z) 



du x 



r 2 {z)^^n(M,z) dz dM 



(31) 



The integration over the spatial correlation function may be performed by substituting the cluster separation R (equation[2£ 
in the dark matter correlation equation 125H . which is in turn substituted in the integral over u in equation 1311 1: 



£ k^T^dk 



sinfc^f ■ u) 2 + (d A (z) ■ 6)2 



du, 



(32) 



N £•« ) +(d A (z) 
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Figure 10. The A2PCF calculated for the Gaussian and x\ models for cases (a) and (b). The left-hand panel describes the correlation 
levels, whereas the right-hand panel depicts the ratios between the x\ an d Gaussian models. 



where the last expression was obtained by changing the order of integration. The integral over u can be performed using the 
following change of variables: 



_ \/(j§ ■u) 2 + {d A {z)-ey 

X ~ d A (z) ■ 6 

which results in an integral that can be put in a closed form 
2 f°° sm[kxd A {z)9]dx _ tt 



k(dr I dz) J j 



-J [kd A (z)6], 



(33) 



(34) 



VPr-T k(dr/dz) 

where x was defined above and Jo is a Bessel function of the first kind and order zero. Expression l|32[l then finally becomes 

„ , , , , ; / k n+1 T 2 (k)J [kd A (z)6]dk. 
2TT(dr/dz) J 

This integral can be easily evaluated numerically. 



(35) 



4.2 Results 

The A2PCFs in cases (a) and (b) in the Gaussian and Xi models are plotted in Fig. 1101 Their shapes are dictated by the 
functional form of the spatial correlation function, which is dominated by the Bessel function of the first kind and order zero 
appearing in equation 13511 . Referring first to the differences between the predictions of cases (a) and (b), it can be seen that in 
both the Gaussian and xi models the correlation levels are higher in case (b). The reason for this is easily explainable noting 
that case (b) is biased towards higher-mass clusters owing to its lack of temperature evolution; as clusters do not become 
hotter with increasing redshift as in case (a), only very massive clusters generate detectable S-Z signals. These higher-mass 
clusters are more strongly correlated, as indicated by the the bias term, equation by virtue of its specific dependence on 
the mass variance <r(M), a monotonically decreasing function of mass M. Thus, the bias factor is a monotonically increasing 
function of M, giving rise to stronger correlations among high-mass cluster populations. 

In order to explain the differences between the predictions of the Xi and Gaussian models in either case (a) or (b), it 
will first be useful to give a qualitative explanation of the shape of the A2PCF curve. We may think of an effective maximum 
physical separation for which (positive) correlation between two clusters is still possible. At low angular separations the scale 
corresponding to this limiting distance can be subtended at relatively high redshifts, so contribution to the A2PCF at this 
angle is possible from low as well as high redshifts. With increasing angular separation the corresponding scale can only 
be subtended at lower redshifts, so that the main contribution to the correlation is generated locally. Since the Xi model 
yields higher cluster populations (and particularly so of high-mass clusters) at high redshifts, it is natural to expect higher 
correlation levels than in the corresponding Gaussian case, where the cluster population is limited to lower redshifts, and 
also relatively limited in massive clusters. As can be seen in the right-hand panel of Fig. 1101 at angular separations of 6 > 3' 
and 6 > 10.5' in cases (a) and (b), respectively, the correlation levels in the Xi model fall below those of the Gaussian 
model. Since the redshift range available for contributing to the angular correlation function reduces with increasing angular 
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separation, this trend is indicative of a diminished contribution of relatively low-redshift clusters in the Xi model. In fact, this 
is indeed the case, and can be easily realized by inspection of the upper and lower- right hand panels of Figs. and (J7J, where 
it has been shown that a substantial contribution to the power spectrum levels at high multipoles and to the cumulative 
number counts originates in clusters lying at high redshifts (z > 1). As explained above, at such high redshifts clusters are 
not strongly correlated at high angular separations by virtue of their mutual large physical distances. On the other hand, in 
the Gaussian model the corresponding contribution comes almost exclusively from low-redshift clusters, which may manifest 
relatively strong correlations even at higher angular separations. Clearly, this effect is also responsible for the steeper shape 
of the A2PCF curve in the non-Gaussian models. 



5 DISCUSSION AND CONCLUSION 

The main objective of this study has been to compare predictions of S-Z observables for two mass functions which differ mainly 
at the high mass end. A Xi distributed density fluctuation field is characterized by a longer tail of high density fluctuations 
with respect to the Gaussian case, and is capable of generating more high-mass clusters at higher redshifts, owing to the 
denser population of high density fluctuations (or peaks) which can collapse at earlier times. S-Z power spectrum, number 
counts, and 2-point angular correlation function all prove to yield significantly different results if calculated using a mass 
function derived from a Gaussian or xi primordial density field. We have demonstrated that levels of the power spectrum 
in the latter model are higher by about an order of magnitude than in the former model, depending on the IC temperature 
scaling with the cluster mass. This was explicitly shown to take effect in light of the larger contribution of high-mass clusters 
and high redshift clusters to the total power in the non-Gaussian model. The larger contribution from distant clusters is also 
responsible for the higher multipoles at peak power. 

We have shown that S-Z power spectrum levels are contributed by different cluster populations in the Gaussian and x\ 
models; in the former model, most of the power originates in clusters lying at < z < 1, with the exception of the higher 
multipoles for which a non-negligible amount of power comes from I < z < 2. In contrast, in the latter model more power is 
contributed from higher redshift clusters, beginning at multipole ~ 3000. On the other hand, most of the power is contributed 
by clusters lying in the mass range 10 14 < M < 10 15 Mq. This holds for both models in the interesting multipole ranges. 

The power spectrum is quite sensitive also to cluster internal properties, as is demonstrated in part through the mod- 
ification of the cluster temperature scaling with mass, and by the dependence on the cluster core radius. These reflect the 
uncertainty of the actual scalings due to observational scatter, lack of high-redshift cluster data, and insufficient theoretical 
information. Plausible scatter of 10 — 20% in these scalings results in significant differences in levels of power, albeit to a 
considerably lower degree than those produced by replacing the Gaussian mass function with the Xi mass function. Number 
counts are also significantly higher in the xi model (up to an order of magnitude) , reflecting the higher population of high- 
redshift, massive clusters, capable of generating sufficient S-Z flux to be detected by upcoming experiments. These results 
obviously depend on the flux detection limit of the experiment, and are likely to become more pronounced with increasing 
limiting flux, owing to the high-mass cluster bias. 

The 2-point angular correlation function of S-Z clusters is potentially an additional observable for which the Gaussian and 
Xi models may result in substantially different predictions. Depending on the temperature scaling with redshift, it has been 
shown that the correlation level at the lowest angular separations (~ 1') is higher by a factor ~ 2 and ~ 4 in the xi model. On 
the other hand, at larger angular separations correlation levels are higher in the Gaussian model. This characteristic behaviour 
leads to a steeper functional form in the xi model. Although not explicitly demonstrated in this study, the flux limit of the 
experiment is, of course, a major factor in determining the correlation levels, due to the fact that higher flux limits are biased 
towards more massive clusters, whose bias parameters are larger. 

An additional relevant property which may be manifested differently in Gaussian and non-Gaussian density fields is cluster 
(or halo) formation time. Clusters are likely to form earlier in (positively skewed) non-Gaussian fields owing to the higher 
probability of finding an overdensity exceeding the critical density for collapse with respect to the corresponding probability 
in a Gaussian density field. This may bear consequences for, e.g., the concentration parameter of haloes, a key parameter of 
the NFW dark matter profile (Navarro, Frenk & White 1995). The concentration parameter is expected to assume higher 
values with earlier formation time by virtue of the higher background density. In this context the recent result from a lensing 
analysis of the cluster AI689 (Broadhurst et al. 2005) is of particular interest. A very steep mass profile was deduced in this 
analysis, one that is fit by an extremely high concentration parameter of c = 13.7lj' 4 . This value seems to be at a considerable 
variance with theoretical predictions. Observations of two other clusters yield c ~ 12 for MS 2137-23 (Gavazzi et al. 2003), 
and c ~ 22 and c ~ 4 for the main and secondary clumps in the cluster CI 0024 (Kneib et al. 2003). Although perhaps not 
very meaningful from a statistical point of view, these results may be indicative of earlier formation times, which are favoured 
in positively-skewed non-Gaussian fields. 

The temperature of the IC gas of virialized clusters may also be affected by the cluster formation time. In conventional 
temperature-mass relations it is the observed redshift of a cluster that determines its temperature. In this recent formation 
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Figure 11. Ratios of effective temperatures are shown. The average temperature predicted in the x| model is divided by the corresponding 
quantity in the Gaussian model for / = 0.5 and / = 0.75. Also shown is the ratio between the effective temperature in the non-Gaussian 
model and the temperature predicted by the temperature-mass relation provided in equation (13). 



approach, the relevant redshift is the observation redshift. Clearly, this picture does not address the merging process that 
gives rise to the clusters we observe today. In the merger picture the formation redshift, Zf, may differ significantly from the 
cluster observed redshift, z obs , and the predicted temperatures will vary accordingly. 

The temperature and concentration parameter dependence on formation time are likely to influence S-Z observables; 
higher IC gas temperatures (stemming from earlier formation times) may give rise to an intensification of the effect, whereas, 
if one assumes that IC gas approximately traces the dark matter distribution, higher concentration parameters may bring 
about transfer of power to higher multipoles, in light of the resulting steeper gas density profile. 

We have conducted a preliminary study of the probability distribution function (PDF) of formation times within the 
framework of Gaussian and xl density fields. This initial work is based on the excursion set formalism and its derivatives 
(Bond et al. 1991, Lacey & Cole 1993). We note that the definition of the formation time used here is somewhat arbitrary; 
clusters were assumed to have formed once they assembled a fraction of either / = 0.5 or / = 0.75 of their current masses. We 
have found that on most relevant mass scales clusters form at higher redshifts in the xl model. Defining "effective" cluster 
temperatures and concentration parameters by averaging these two quantities at specified formation redshifts over the PDF of 
formation times, we learn that temperatures and concentration parameters are on average higher in the xl model for masses 
10 13 - 10 15 M Q . At higher masses the differences are practically negligible. Also, with increasing observation redshift, the 
differences become gradually less pronounced. This is related to the fact that the time range available for formation is reduced 
with increasing observation time, since the cluster must have formed at a redshift exceeding the observation redshift. 

As an example of our preliminary results we include Fig. il l H , which compares the effective IC gas temperature (as defined 
above) in the xl and Gaussian models as a function of cluster mass, for / = 0.75 and / = 0.5. Also shown are the respective 
temperature ratios between the xl predictions and those given by equation 1131 . Clearly, the indicated average temperatures 
are closely related to the redshift range at which the formation time PDF attains the highest levels. For example, the PDFs of 
the lowest mass scale (1O 13 M0) peak at relatively high redshifts, and therefore the weight of the contribution from the (1 + z) 
factor is rather large. This results in higher average temperatures. The effect is most pronounced when comparing the xl 
model with the non-averaged temperature of equation l|13|l . since for the latter case all clusters are assumed to have formed 
at Zobs = 0. Obviously, / = 0.5 implies earlier formation times and higher average temperatures for a given (low) cluster 
mass. This holds also for the xl model, which involves earlier formation times for this mass range. The average temperature 
predictions in the xi an d Gaussian models show that clusters of 10 13 Mq are hotter by factors ~ 1.4 and ~ 1.27 for / = 0.5 
and / = 0.75, respectively. For M = 10 14 Mq the corresponding factors are ~ 1.23 and ~ 1.12, while for M — 10 15 Mq in 
both cases the factors are ~ 1. 

We may therefore conclude as follows: a mass function based on a Xm density fluctuation field gives rise to higher 
populations of massive clusters at relatively high redshifts with respect to a mass function derived on the basis of Gaussian 
density fluctuation fields. This generates higher levels of S-Z power spectrum and number counts by up to an order of 
magnitude. The 2-point angular correlation function of S-Z clusters is also susceptible to the mass function model and predicts 
higher correlations at low angular separations in the non-Gaussian model. At higher angular separations the correlation levels 
in a Gaussian field dominate. Cluster formation times are generally earlier in the non-Gaussian model, and may also give 
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rise to higher temperatures and concentration parameters of low-mass clusters observed at relatively low redshifts. A more 
rigorous exposition of the topic of formation times, their relevance to S-Z observables, and further related ramifications will 
be discussed in a future paper. 
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